Langevin description of nonequilibrium quantum fields 
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We consider the non-equilibrium dynamics of a real quantum scalar field. We show the formal 
equivalence of the exact evolution equations for the statistical and spectral two-point functions with 
a fictitious Langevin process and examine the conditions under which a local Markovian dynamics is 
a valid approximation. In quantum field theory, the memory kernel and the noise correlator typically 
exhibit long time power laws and are thus highly non-local, thereby questioning the possibility of a 
local description. We show that despite this fact, there is a finite time range during which a local 
description is accurate. This requires the theory to be (effectively) weakly coupled. We illustrate 
the use of such a local description for studies of decoherence and entropy production in quantum 
field theory. 



I. INTRODUCTION 

Understanding the dynamics of time evolving quantum 
systems is a key issue in many topical areas of physics, 
from early-universe cosmology to high-energy nuclear col- 
lisions, condensed matter physics or ultracold atomic 
gases. One important line of investigation in nonequilib- 
rium field theory concerns the issue of first principle cal- 
culations of far-from-equilibrium (quantum) dynamics. 
Major breakthroughs have been achieved in recent years, 
in particular with the use of two-particle- irreducible func- 
tional techniques |l|, [2[ , and the field is under active de- 
velopment, thanks to the advent of ever faster computers 
and of new ideas 

An important byproduct of such investigations is that 
it allows one to bridge the gap between relatively simple 
situations where direct calculations from the basic equa- 
tions of quantum field theory (QFT) can be done and 
more intricate cases, such as, for instance, situations with 
expanding and/or inhomogeneous backgrounds or 
situations involving many types of fields and interactions 
such as lepto/baryo-genesis scenarios , where such 

calculations are difficult and where one often has to rely 
on effective descriptions. The most popular such effec- 
tive approaches are kinetic descriptions [1, HI, fT3 - [20| | , or 
effective Langevin descriptions (2l| - [28j . A key issue in 
this context concerns the domain of validity of such ap- 
proaches. 

The derivation of kinetic - e.g. Boltzmann - equations 
from nonequilibrium Schwinger-Dyson, or Kadanoff- 
Baym equations have been extensively discussed since 
the early works of Kadanoff and Baym [|, [H-El [H [H- 
l33j . This typically relies on a gradient expansion, which 
assumes a clear separation between the scale of the 
nonequilibrium dynamics at hand and that of the elemen- 
tary excitations and processes responsible for it. The va- 
lidity of the gradient expansion and of the corresponding 
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transport equations has been studied against direct solu- 
tions of the underlying nonequilibrium Schwinger-Dyson 
equations in cases where the latter can be done |34l439{ . 

Effective Langevin equations and their derivation from 
basic QFT have also been much studied, mainly along 
the lines of the pioneering work of Feynman and Ver- 
non, based on the so-called influence functional formal- 
ism [23|, [25, 40 44]. The latter provides a systematic 
way of integrating out irrelevant degrees of freedom in a 
path integral language. One typically obtain non-local - 
i.e. non-Markovian - Langevin-like equations with vari- 
ous additive and/or multiplicative random noises |4l[|45l |. 
Noise correlators and memory kernels can in principle be 
computed by means of standard QFT techniques. 

The most popular use of Langevin descriptions involves 
further drastic simplifications, namely the assumption of 
effective local damping and Gaussian white noise, result- 
ing in a Markovian dynamics (2lT[28j . This assumes a 
clear separation of scales between the so-called "memory 
time" , usually identified with the short time scale of ir- 
relevant degrees of freedom, and the "relaxation time", 
which characterizes the dynamics of the relevant degrees 
of freedom [46l . . Existing comparisons between given 
non-Markovian L ang evin dynamics and their Markovian 
versions, see e.g. (4a-l5lj indicate that, indeed, the better 
the separation of scales, i.e. the weaker the local damp- 
ing, the better the Markovian approximation. 

As mentioned above, memory kernels can in princi- 
ple be computed in a given theory from, say, loop di- 
agrams. Existing calculations in simple theories with 
scalar, fermionic and/or gauge fields show that memory 
kernels and noise correlators typically decay as power law 
in time, thereby questioning the existence of a memory 
time and more generally the possibility of a local, Marko- 
vian description [H, [Hj]. Of course, there are situations 
where one expects a local Markovian Langevin descrip- 
tion to be questionable, such as in presence of massless 
(e.g. Goldstone or gauge) excitations and/or at zero tem- 
perature (53l . Ilij . But the absence of a clear memory time 
scale mentioned here also occur in more standard situa- 
tions with no infrared issues and at high temperatures, 
where one would expect a Brownian like motion to be a 
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good description. 

We address this issue in the present paper. We first re- 
call the basic evolution equations for nonequilibrium two- 
point Green's functions, which involve non-local memory 
kernels (self- energies), and show their formal equivalence 
with a fictitious Langevin process. We discuss the basic 
conditions under which the memory integral and noise 
kernel can be replaced by local - mass and damping - 
terms and argue that these conditions do not require the 
existence of a local limit of the kernels themselves. In 
order to study the validity of a local Markovian limit 
of the equivalent fictitious Langevin process, we focus 
on a simple situation where an out-of-equilibrium test 
field is coupled to a thermal bath with negligible back- 
reaction. The basic QFT equations can then be exactly 
solved in terms of the equilibrium spectral function of the 
test field. We find that there exist a characteristic time 
below which the actual dynamics is indeed Markovian 
but after which memory effects cannot be neglected, as 
expected on rather general grounds (5^-[H3 • We investi- 
gate this time scale in detail in a simple model with cubic 
interactions between the test and the bath fields. Finally, 
we discuss the consistency of the local limit directly in 
real time, at the level of the equations of motion. 

As an illustration, we end this paper with an appli- 
cation of such local description to the physics of deco- 
herence and entropy production in QFT in the so-called 
incomplete description picture recently advocated in (Hs" 



II. THE STRATEGY 
A. General setting 

We consider a generic Z2-symmetric scalar field the- 
ory in the symmetric phase. The nonequilibrium n- 
point functions can be conveniently described by means 
of time-ordered products of field operators on a closed 
contour C in time (66l The two-point function 

G(x,y) — (Tc(f(x)(p(y)) encodes both the statistical 
and spectral correlators F(x,y) = ^({ip(x),cp(y)}) and 
p(x,y) = i{[<p(x),<p(y)]) as 

G(x, y) = F(x, y) - ^sign c (a; - y°)p(x, y) . (1) 

Here, the brackets (...) denote an average with respect to 
a given density matrix describing the (quantum) state of 
interest. The common practice is to specify the (out-of- 
equilibrium) initial conditions at a given finite time t = 
[if, HI ■ In the present work, we choose instead to prepare 
the system in a Gaussian thermal state at t — > — oo and 
to kick it away from equilibrium at t = by means of 
an appropriate external source. Both procedures are in 
principle equivalent. The source approach is an efficient 
device to let the system develops its own non-Gaussian 
correlations [68:1, 69]. The latter are of particular impor- 
tance e.g. when it comes to discussing renormalization - 



the ultraviolet modes have to be in the correct non-trivial 
vacuum state [69|, [zO] ■ We shall not be concerned with 
these issues in this paper and adopting the source ap- 
proach is only a matter of technical convenience for later 
calculations. Still we present it in some detail because of 
its potential interest also for numerical calculations. 

We consider a Gaussian nonequilibrium disturbance 
described by a bilocal source K in the action [71| : 

S K [<P] = S[<p] + \( dkd*u<p{x)K(x,v)<p{y) , (2) 
J c 

where S is the classical action of the system. The 
Schwinger-Dyson equation for the two-point function on 
the contour C reads: G _1 = G^ 1 — £ + K, where 

iGQ 1 (x,y) — — (□ + M 2 )8^\x — y) is the free propa- 
gator for the field of mass M and Sfajj/) its self-energy. 
Using the standard decomposition [lj, |2j 

E(x, y) = S (a;)4 4) ( x ~ v) 

+ Z F (x,y)- ^sign c (x -y°)Z p (x,y). (3) 

and assuming Gaussian conditions at t — > — oo, the 
Schwinger-Dyson equation translates into the following 
nonlinear coupled integro-differential equations for the 
statistical and spectral propagators 

r° 

[n x + M 2 {x)]F{x,y) = - / d 4 zE p (i,z)F(z,y) 

J — OO 
n 

d 4 z^(x,z)p(z,y),(4) 

r° 

[a x + M 2 (x)]p(x,y) = - d 4 zZ p (x,z)p(z,y), (5) 

JyO 

where d 4 z = J Q b dz° Jdz, M 2 (x) = M 2 + E (x) and 
V$(x,v) = H F (x,y)-K(x,y). (6) 

For later purposes it proves convenient to rewrite Eqs. 
(HI)-© by introducing the retarded and advanced prop- 
agators G R (x,y) = 9(x a - y a )p(x,y) and G A (x,y) = 
Gn(y, x) — — 6(y° — x°)p(x, y) and similarly for the self- 
energies. One gets, from Eq. ([5]), 

G^( Xi y)= p x + M 2 (x)]sW(x-y) + Z R (x,y) (7) 

and similarly for GA{x,y) = Gn(y 7 x). Eqs. (HJ)-© can 
be written 72] 



F = -G R *(T, F -K)*G A , 
P = —Gr * E p * Ga , 



(8) 
(9) 



with (f*g)(x,y) = Jd 4 zf(x,z)g(z,y). Eqs. 0-© ac- 
tually provide an explicit solution if the self-energies are 
known. In general, however, the latter are non-linear 
functions of the two-point functions themselves. 
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The source K serves to prepare out-of-equilibrium con- 
ditions near x° = y° = 0. For instance, we shall consider 
an instantaneous kick at x° = y° = 0: 

K(x, y) = A(x, y)5(x°)S(y°) + B(x, y)S'(x°)S' (y°) 

+ C(x, y) [6'(x°)6(y°) + S(x°)S'(y )] . (10) 

Before the kick, the system has had an infinite amount of 
time to reach an equilibrium state with parameters (e.g. 
temperature) determined by the conditions (e.g. energy 
density) at time t —¥ — oo: 



F ( x >v)\xo,yO <0 = F eq (x-y). 
Right after the kick, we write: 

F(x,y)\ x o ty0>0 = F ncq (x,y). 



(11) 



(12) 



The functions A, B and C in (JTDJ) are related to ini- 
tial conditions for the nonequilibrium statistical function 
^neq [51 • From Eq. ©, we obtain 

F(x, y) = - [G R * E F * G A ] (x, y) + 6(x°)9(y°)F K (x, y), 

(13) 

with 



F, 



K 



(x, y)= {p(x°, 0)-A-p(0, y°)+d u p(x°,u)-B-d v p(v, y°) 
-d u p(x°, u)-C-p(0, y°) - p(x°,0)-C-d v p(v, y°)\ 

J u—v—0 

(14) 

where the dot denote a convolution with respect to spa- 
tial variables only. Using d x op(x,y)\ y o =x o = c)( 3 )(x — y) 
we obtain the relations: 

B(x,y) = F K {x,y)\ x0=y0=Q , 

C(x, y) = d x oF K (x, y)\ x0=y0=0 , (15) 

M x >y) = d x od y oF K (x,y)\ x0=y0=Q . 

Somewhat more intuitive expressions of the sources can 
be obtained by defining 



and 



<(x,y) 
<(x,y) 
^ c ° q (x,y) 

^„ cq (x,y) 

-Rneq( X >y) 

^„° cq (x,y) 



F[X, yJla-o^o^o- ' 

d x oF(x,y)\ x0 y0 ^ - , (16) 
d x od y oF(x,y)\ x0y0 ^ _ 

F(x, y)\ x o jy o^ + , 

d x oF(x, y)\ (17) 
d x od y oF(x, y) | x o o_> 0+ ■ 



Assuming that the function Gr * E F * G a is continuous 
at x° = y Q = 0, one has 



B(x,y)=F n %(x,y)-F ° q (x,y), 
C(x,y)=i?° cq (x,y)-i?° q (x,y), 
A(x,y)=X° cq (x,y)-X o q (x,y). 



(18) 



B. Equivalent Langevin process 

The evolution equations (0-® are formally equivalent 
to a fictitious Langevin process <p* driven by a random 
noise £ [z3 : 

[D x + M 2 (x)] ip 6 {x)+ J d 4 z£ R (x,z)<pt{z) = £(x)+E(x), 

(19) 

where 2 is a fluctuating disturbance source which serves 
to send the system away from equilibrium at a given finite 
time, in the same spirit as the source K above. The 
random variables £ and E are uncorrelated. 

Using 0, Eq. can be rewritten as G^ 1 *^ = 
whose solution is formally given by: 

/(:r) = <A)(z) + / d 4 z G R (x, z) [£(*) + E{z)] , (20) 



where tpo is the general solution of the homogeneous 
equation G^ 1 * tpo = 0. 

We shall denote the average over realizations of the 
noise £ with an overline and the average over the source 
5 with double brackets. The original quantum average 
is identified with (...) = ((. . .)). Setting ((E{x))) = 
and £(x) = and choosing initial conditions such that 

<Po{x) — guarantees that (y(x)) = ((ip^(x))) = 0. To 
obtain the evolution equation for the two-point function 



we write (fT9")l as 



(21) 



d 4 zG- 1 (x,z)^(z)^(y) = [ax)+E(x)]^(y). 

. (22) 

Inserting Eq. (1201) on the right-hand-side and identifying 



Z F (x,y) = -Z(x)Z(y) and K(x, y) = «H(x)S(y)» , 

(23) 

one recovers Eq. flSJ). 

The equivalence is complete provided the sources are 
chosen in an appropriate way. For the case of an instan- 
taneous kick, 

H(a;) = r(x)<5(x°) + n(x)<5'(x°) , (24) 

one gets 

B(x,y) = ((n(x)n(y)», 

C(x,y) = i«T(x)n(y) +n(x)T(y)}) , (25) 

A(x,y) = <(T(x)T(y)>), 

with ((T(x)n(y) -n(x)T(y)» =0. 

It is worth emphasizing that, in general, the noise cor- 
relator Ef and the memory kernels E p depend nonlin- 
early on the correlators F and Gr and are thus nonlin- 
ear functions of the field yfi and the noise £ and their 
derivatives. In particular, this encodes all sorts of addi- 
tive and/or multiplicative noises as well as linear and/or 
nonlinear damping terms that one encounters in usual 
(mostly perturbative) derivations of Langevin equations 
e.g. based on the influence functional technique |25l. [26|. 
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C. Approximation strategy 

We now examine sufficient conditions under which the 
memory integral in the exact nonlocal Langevin equation 
(Tl9|) can be approximately described by a local damping 
term. We shall see that there is no need a priori for as- 
suming a finite memory time. For simplicity, we shall 
specialize to spatially homogeneous and isotropic situa- 
tions. In particular the source K and, consequently, all 
other two-point functions in the problem can be spatially 
Fourier transformed: 



K(x 



d 3 p 



\ f d > 



-^K p (x°,y°) (26) 



and similarly for all other two-point functions. The rele- 
vant equations read 



[d 2 +^ 2 (i)] F p {t,t') + / dumt,u)F p (u,t') 

J — OO 

X) 

duY,^ K {t,u)Gj{u,t') (27) 



and 



[df + cj 2 p (t)] G*(t,t')+ / duEg(t,u)G£(tt,f) = 6(t-t'), 

J — OO 

(28) 

where w p (t) = y/p* + M 2 (t) and T,^ K = - K p . 

In generic theories, the memory kernel t') is typi- 
cally a decreasing function of t — t' . For instance, writing 
the interaction term of the field degree of freedom of in- 
terest ip p (t) as 5i„t = f t p ip p (t)Jp(t), where the current 
J p involves other degrees of freedom, possibly including 
tpp'^p, one has, at lowest order in the interaction, 



mt,t')<x{[j p (t),j p (t')}) 



(29) 



Such current-current correlator typically exhibits rapid - 
e.g. power law [H, [H, [75[ - decay at large t — t' and the 
memory integrals in Eqs. (|27|l and (|2"B"|) are dominated 
by their upper bound. 

As for the functions F p and G p , they are generically 
characterized by various - possibly time-dependent - 
time scales: an oscillation frequency e p (i), the damping 
scale in the relative time t — t' and the typical scale of the 
nonequilibrium evolution of equal-time correlators - the 
last two might actually be power laws. The most common 
situation - e.g. for weakly coupled (effective) theories - 
is that the former is the shortest time scale in the prob- 
lem. We shall assume a clear separation of scales, where 
the oscillation frequency is always much shorter than the 
damping or nonequilibrium scales. 

Now, if the range of integration which dominates 
the memory integrals is short compared to the above- 
mentioned damping and nonequilibrium scales, one can 
neglect the latter under these integrals and replace the 
full functions F p and G p by their short-time oscillating 



part around the upper bound of the integral, e.g. |76j |: 
F p (u,t') w F p (t,t') cos ep(t)(u-t) 

_w^n sinepmu - t ), (so 
e p\ i ) 

with some relative error controlled by the ratios of either 
the damping or nonequilibrium scales over the oscilla- 
tion frequency. The memory integral in (127[) can thus be 
approximated as 



OME p p (t,u)Fp(u,t') « 5<? p {t)Fp{t,t'Y-2lp{t)d t F p {t,t'), 

(31) 



with 



and 



5e 2 p (t)= Reflet; e p (t)) 
Im£*(t;e p {t)) 



7p(*) 



2e p (t) 



(32) 



(33) 



where we used the definition E^(i, t') = 6{t - t')T,P(t, t') 
and introduced the mixed time-frequency representation 



E*(t,t') 



2^ 



e -M*-0£fl(i. w ). 



(34) 



Similar manipulations hold for the memory integral in 
Eq. (EH). 

If the above conditions are met, the equation for F p 
and G p thus read 

[d* + 2 7p (t)d t + el(t)]F p (t,t') = 



+ 00 



AuY>^ K {t 1 u)G p x {u,t') (35) 



and 



[d 2 t + 2 lp (t)d t + e 2 p (t)] G«< A (t,t>) = 5(t - , (36) 

where we defined e 2 (t) = ui p (t) + Se 2 (t). Using similar 
manipulations as in the previous subsection, it is easy 
to check that these equations also have an equivalent 
Langevin process which reads: 

[d 2 + 2 lp {t)d t + e 2 p {t)} 4(t) = £ p (t) + E p (t) . (37) 

Note that at this stage, the memory kernel has been re- 
placed by a local damping term but the noise correlator 
can still be non-trivial, i.e. non-Markovian. Note also 
that, here, the damping rate and the noise correlator are 
still, in general, complicated non-linear functions of the 
propagators. 



III. EXACTLY SOLVABLE EXAMPLE: NEAR 
STATIONARY SYSTEM 

We now want to investigate in more details the condi- 
tion under which the above local approximation is valid. 
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In this section we specialize to a simple situation where 
the out-of-equilibrium degrees of freedom weakly interact 
with a stationary - e.g. thermal - background with neg- 
ligible backreaction. In this case, the full nonequilibrium 
dynamics can be exactly solved [H, [H, [44j and one can 
check whether the exact solution can be described by a 
Langevin dynamics with local damping. 



A. Exact solution 

The assumption of negligible backreaction means that 
the self-energies are determined by the station- 

ary background and are thus time-translation invariant: 
Ep' p (t, t') = ^'P(t-t'). It follows that the spectral func- 
tion is also time-translation invariant: p p (t, t') — p p (t—t') 
[l~8j . Indeed, using the time-translation invariance of ££, 
one easily checks that the equation of motion for p p (t, t 1 ), 
see e.g. (|28|) . only explicitly involves the time differ- 
ence t — t' and thus possesses time-translation invariant 
solutions. Moreover, since the equal-time commutation 
relations, which determine the initial conditions for p p , 
are preserved under time-evolution, the only possible so- 
lution has this symmetry. The equation of motion for 
p p (t — t') reads: 

[d 2 t + uP p ] Pp (t) + [drllftt- t) Pp (t) = , (38) 
Jo 

with p(0) = and p(Q) = 1. It follows that the functions 
Gp' A and, in turn, F cq = — Gr * T,p * Ga are also time- 
translation invariant. 

Introducing the frequency representation 



implies 



P P (t) 



2^ 



(39) 



and similarly for all other time-translation invariant func- 
tions Gp' A and Ep' p ' fl ' yl , one readily obtain 



G» = 



as well as 



-(w + iO+) 2 + w2 + £fi( w + iO+) 

f;^) = -|g«h| 2 sJh, 

p p (u) = -\G«(uj)\ 2 £p'(u). 



(40) 

(41) 
(42) 



It is useful to recall that p p (uj) = 2ilmG^(uj) and sim- 
ilarly for the self-energy components [73]. Notice the 
detailed balance relation 



(43) 



characteristic of stationary systems. In the following 
we assume a thermal background to fix the ideas al- 
though the argument does not depend on this assump- 
tion. In thermal equilibrium at temperature T = 1//3, 
the fluctuation-dissipation relation 



(44) 



(45) 



with n(x) — l/(exp(x) — 1) and thus there is only one 
independent function, e.g. p p (uj). 

For an instantaneous kick at t = t' = as in (fTUj) : 

K p {t,t') = ApS(t)S(t')+B p S'(t)5'(t') 

+ C p [8'(t)6(t')+S(t)S'(t')] , (46) 

the complete nonequilibrium solution for F p reads, see 

Eqs. (ni- dm, 



F p (t, t') = F^(t - + 6(t)9(t')Fp K (t, t') , (47) 



with 



where 



F*(f,i') = A pP p(t)pp(t') + Bppp(t)pp(t') 

+ Cp[pp(t)pp(t')+p p (t)pp(t')] , (48) 



Bp - (*,*') | t=t ,= , 
C p = d t F p «(t,t>)\ M 



(49) 



d t d t ,F«(t,t' 



Thus we see that, in the simple situation considered 
here, the nonequilibrium dynamics is fully described by 
the second term on the right-hand-side of Eq. (|4"7} : The 
momentum modes of interest get kicked away from equi- 
librium at t = and the relaxation toward equilibrium, 
described by Eq. (|4"8"|) . is completely encoded in the (equi- 
librium) spectral function p p (t). 



B. Breit-Wigner approximation 

In general, the spectral function p p (uj) may exhibits 
single-particle poles as well as branch cuts in the complex 
frequency plane, correspondin g to multiparticle processes 
or Landau damping, see e.g. |52l. l53l. . For instance, 
near the poles, it takes the usual Breit-Wigner form [79| 



p p (u) -> e p ) 



*1p 



€p [U> — c p 

where we assume 7 P /e p <C 1. Here, 



<r,) 2 +l 2 p 



(50) 



(51) 



and the quasiparticle energy e p and width 7 P are deter- 
mined from 



and 



el=uj 2 p+Re^(ep) 
Imfff(e p ) 



7 P = -Z p - 



2e r , 



(52) 



(53) 
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Poles lead to exponential damping of the Fourier trans- 
form p p (t) at large time, whereas branch cuts result in 
power law behavior [EM Hit- The large-time behavior is, 
typically, 



p P (t) 



sin^ e _ 7p|t| 



(54) 



where the last term denotes the power law contribution. 
Using (|44|) one can obtain the corresponding large-time 
behavior of the equilibrium statistical function F° q (i— t'): 

Fp^t) « K P (cose p t + 2e sme p \t\^ + FP° w (t) , 

(55) 

with 



n (Pep) 



(56) 



The power law contributions to both p p (t) and F^(t) 
are due to features of the self-energies £^ p (u;) and their 
amplitudes are thus governed by the strength of the inter- 
action with the thermal bath, i.e. by some power of the 
relevant coupling constant. This is to be contrasted with 
the pole contributions, where the interaction strength es- 
sentially controls the decay rate j p but only gives correc- 
tions to the amplitude. Thus, for weak coupling, there 
exist a range of time during which the exponentially de- 
caying contributions dominate over the power laws. Only 
at very late time is the dynamics described by the latter. 

Writing schematically the power law contributions as 
~ A/(jJ,t) u , with A an amplitude proportional to some 
positive power of the relevant coupling constant or, equiv- 
alently, of the damping rate j p , p some scale of the prob- 
lem at hand and v a given exponent [80], one finds that 
the power law contributions arc of relative order jp/e p 
for times 



7 P 



*<ln Z- 



4A 



(57) 



In this time range, the nonequilibrium relaxation 
also essentially exponential: 



is 



F p K {t,t')^F°{t,t')e-^ t+t '^ 



(58) 



with 



F°(t, t') = a+cos e p (t-t')+a~ cos e p (t+t')+a° sin e p (t+t'), 

(59) 

where the constant and a° are related to initial con- 
ditions dm} 



7 2 I 
a ± = fjL Bv ± 



A p - 2-fpC p + j 2 B p 



Thus the complete nonequilibrium solution takes a re- 
markably simple form in the time range (|57| 81]. In the 
weak coupling limit both jp/pp and e p A p are small and 
the latter can be relatively long. In that case, the de- 
viation from equilibrium (|58[) as well as the equilibrium 
unequal-time correlators (|54|) and (|55|) have essentially 
completely decayed when the power law contributions 
start becoming important and most of the relevant dy- 
namics is thus well described by the exponential - pole - 
contributions. As we now discuss, the latter can be de- 
scribed by an equivalent Markovian Langevin dynamics. 



C. Markovian Langevin dynamics 

In the time range (|57j> . the dynamics is dominated by 
the poles of the retarded Green's function (|4"0"|) . Those 
which are closest to the real axis give the dominant late- 
time contribution. Assuming that the dynamics is driven 
by these poles and that the Breit-Wigner approximation 
(I5TJ1) is justified, one can make the replacement G p (uj) —> 
,» with H 



-( W + i 7 p) a + na ! 



(61) 



where £l 2 = e 2 — j 2 f=s e 2 . The corresponding spectral 
function reads 



Peff.pM 



■iiZpjpUJ 



{iO 2 - £ 2)2 + 4 7 2 W 2 ' 



(62) 



which reprocuces Eq. (|50p near the poles. The Fourier 
transform reads 



(63) 



The function (|6ip is the frequency space represen- 
tation of the retarded Green's function of the opera- 
tor Z- 1 [df + 2~/pd t + ep] S(t - t'). It follows that the 
nonequilibrium relaxation (|58[) can be described by a 
Langevin process with local damping: 

[8 2 t + 2 lp d t + e 2 p ] $(t) = + Sp(*) (64) 

with appropriate sources S p (t). The equilibrium state is 
characterized by the £ p (t) correlator. A Gaussian white 
noise 



£ v (t)£ p (t>) = Z- 1 a p 6(t-t') 
gives, after simple calculations, 



Z p a p 



(cj 2 



4ry 2 u; 2 



(65) 



(66) 



or, equivalently, 



_ ry2 Cp 1 P B P 

a p - L v ; 



(60) 



, ( cose p t H sine p |i| ) e" 

^lp e P 



-7pI*I 



(67) 
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Choosing 



n„ = 4j p £ p ( n(f3e p ) + - 



(68) 



one reproduces the exponential term in (|55|) . Thus the 
simple Langevin dynamics (|64[) -(|65 p with local damping 
and white (Markovian) noise related by (1551) correctly de- 
scribes the full nonequilibrium dynamics of the previous 
subsection in the time range (|57p . 

In the present case, both the local damping and the 
Markovian nature of the noise follow from the Brcit- 
Wigner approximation and the assumption of a ther- 
mal bath. Indeed the Breit-Wigner ansatz (pjTj) . see 
also (|5U|) . which leads to local damping ~ 7jA, im- 
plicitly assumes that E£(u;) is sufficiently smooth on a 
range \e p — u\ < j p . For the case of a thermal bath, 
the fluctuation-dissipation relation (J45j) guarantees that 
T,p (lj) has the same property. Assuming that the poles 
give the dominant contribution to Fp q (t) amounts to ne- 
glecting the frequency dependence of E^ around the pole 

and replacing E p (w) » E p (e p ) in Eq. (|4T|) . In terms of 
the equivalent Langevin description, this corresponds to 
white noise. We see that, in principle, one may have local 
damping but non-local, colored noise e.g. in the case of 
some non-thermal background. 

Finally, let us recall that the relation (|6"8"|) follows from 
the on-shell fluctuation-dissipation relation (|4~4"| [23| . In- 
deed, the effective Langevin dynamics corresponds to the 
following effective self-energy kernels: 



yR 



p(w) = <5e 2 



2iZl 



7pW. 



with Sep = ReTp(e p ) and 



J cff 



pH = -z p 



S^(e P ). 



(69) 



(70) 



In particular, one has Eg ffp (w) = — 4«Z p 'jpU) and thus, 



-shell, 



(71) 



Eq. (|68j) then directly follows from (|45l taken at ui = e p . 
It generalizes the standard relation a v = 47 P T valid at 
high temperature /3e p < 1 [H |H Equivalently, 
the equilibrium correlators (|62|) and (|67| also satisfy 
the fluctuation-dissipation relation on-shell. Here, the 
fluctuation-dissipation relations (j4"4" |) - (|4"5)) only hold on- 
shell, as a consequence of the assumption that the dy- 
namics is governed by the poles. 

To close this subsection, let us remark that the effective 
memory kernel and noise correlator are completely local: 



Srff, P ( < ) 



-Z p 1 otp5(t) 



and 



(72) 



(73) 



As explained in the introduction, this seems in contradic- 
tion with actual calculations of such functions in various 
models, which typically give power laws in time [53.153$. 
We present an explicit example in the next section. Then, 
we study in Sec. |V]how the apparent contradiction is re- 
solved. 



IV. AN EXPLICIT EXAMPLE 

To illustrate the point , we consider a simple model with 
non-trivial dissipation 0, HH [53[ where the (thermal) 
background is represented by a scalar field x °f mass m 
interacting with the system field ip via a ipx 2 interaction: 

s = - Jd A x^(n + M 2 ) v + 9 -v X 2 

+ i X (D + m 2 ) X + V7x)j. (74) 

We assume weak enough interaction between the two 
fields and strong enough x se ^ interactions such that 
back reaction can be neglected and x remains in ther- 
mal equilibrium at temperature T. We shall compute 
the relevant self-energies in perturbation theory at low- 
est non-trivial order, assuming a free thermal gas for the 
X field. The spectral and statistical components of the 
space-time translation-invariant self-energy of the system 
field Yi[x,y) = E(x — y) read, at one-loop, 



s p(z) = -g 2 F x (x)p x {x) 



E F (x) = - — 



(75) 
(76) 



where F x and p x denote the free statistical and spec- 
tral two-point function of the field \, respectively. Their 
spatial Fourier transforms read (uj v 



F*(t)= n( j 9w p ) + - 



cos ujpt 



vV 



sin ujpt 



(77) 

For simplicity, we consider only the self-energies (|75[1 - 
(1761) at zero momentum in this section. The discussion 
is similar for non- vanishing momentum but is a bit more 
involved due to Landau damping effects which contribute 
extra features in the complex frequency plane. 

Before going on, we mention a peculiarity of the 
present model with cubic interaction vertex: the zero 
mode of the statistical self-energy Ep =0 (t) acquires a 
time-independent contribution and thus does not decay 
to zero at large times. This is most easily seen directly 
in the mixed time-momentum representation: 



sLo(t) 



H 

' 2 



rf 3 



q 



(2tt) 3 



[*?(*)] ! 



f(<)] 5 



The constant contribution arises from the first term 
on the right-hand-side, using ([77)) and 2cos 2 u;,j£ = 
cos 2uj q t + 1 [83j . The actual value of this constant is 



8 



not important for the present discussion and we refer the 
interested reader to Appendix [X] for more details. In 
the following we denote it by E^ =0 (i — > oo) = (non 
constant contributions vanish at large time, see below). 
In the rest of this section, we omit the index p = for 
simplicity. 

The self-energy (|75l) - (|76p read, in momentum- 
frequency space [4J, l78fl 



t p (uj) = -^-e{u? -Am 



Aw 



1 - 



4m 2 



/3oj 



1" 

2 

and (see Appendix [X] for a discussion of the S(uj) term) 



considerations in the previous section, there is a well- 
defined time regime in which the two-point correlators 
F(t) and pit) can be accurately described by a Marko- 
vian Langevin dynamics with local kernels. To make the 
argument more precise in the present model, we now com- 
pute the explicit late-time behavior of these correlators. 
Again this can be done by standard contour integration 
techniques using Eqs. (|4"T|> and (J22J. The calculation 
is actually rather similar as the previous one for self- 
energies with two more pairs of complex conjugate poles 
in the complex frequency plane coming from the retarded 
and advanced propagators Gr(oj) and Ga(oj) — G* r (uj), 
located at ±M - i-y and ±M + ry [si|. with 



E F (w) = -i ( nifiu) + - ) E p (w) + a^Siu) . (80) 



In Eq. (|79|). one recognizes the two particle threshold at 
u> > 2m which implies that the on-shell damping rate 
is non-zero only if M > 2m. The square root factor 
is the standard threshold function for a two body decay 
and is responsible for the branch cut singularity in the 
complex frequency plane, which governs the late time 
power law behavior. Finally the factor n(/3u>) + 1/2 is 
due to Bose enhancement and yields poles at imaginary 
frequencies P<jj n — Aniw, with n € Z. 

The late time behavior of the Fourier transforms 
^•F,p{t) can be obtained by standard contour integration 
techniques. For m ^ 0, the regime Tt 3> 1 is dominated 
by the t wo p article threshold at w = 2m. We obtain, see 
also HIM III, 



E p (t)^ ( x^(t)=-a p 



cos(2m< + 7r/4) 
(mi) 3 / 2 ■ 



v u\ ~ pow m , sin (2nd + tt/4) 

hF{t) « CTp (t) + CTqo = Q F } l-foo, 



(mi) 3 / 2 



with 



o 2 m ( I 



n(2(3m) + - ) a p 



(82) 



(83) 



g 2 m 
16tt 3 / 2 



n (/3m) 



The first line of Eq. is reminiscent of the fluctuation- 
dissipation relation (|45| and the second line uses the iden- 
tity 



n(2x) + i 



, x 1\ 1 
n (x) + - = - 
y } 2 2 



n(x) + - 



As announced, the memory and noise kernels E p (t) 
and (t) exhibit highly non-local power law behaviors 
at large times. Still, as explained from rather general 



it p {M) _ g 2 



AM 



16ttM 



4m 2 
M 2 " 



(86) 

The late-time behavior, Tt ^> 1, is governed by these 
poles and by the branch cuts \ui — 2m\ > along the real 
frequency axis. If the poles are sufficiently far from the 
threshold uj = 2m, the pole and cut contributions are 
well separated. Here we assume M 3> m for simplicity. 
We get, for the spectral function, 



p{t) „ ^m e -7i«i 



pP° w (t) (87) 
and, for the equilibrium statistical one, 



F eq (t) « n (j3M) + 



1\ cos Aft 



e-^+F pow (t) + F a 



M 



In both cases, the first term on the right-hand-side is 
due to the poles whereas the second ones arise from the 
branch cuts near threshold: 



pP™(t) = -\G R (2m)\ 



'(*) 



'it) 



(M 2 - Am 2 ) 2 



(89) 

and similarly for F^(t) with crP ow (t) -> cr£ ow (i). Fi- 
nally the statistical correlator does not decay to zero at 
large t because of the constant <Joo in (|82|) : 



^oo =CToo|Gfl(0)| 



M 4 ' 



(90) 



(84) In Eqs. (|59")) - ([W|) . we used the exact expression 



\Gr{u>)\- 



M 2 + ReE fl (w) 



u 2 ] 2 + jS», (91) 



as well as the fact that at the two-particle threshold 



E p (2m) = 0, sec Eq. (|79|1 . 

We are now in a position to obtain a more precise es- 
timate of the time range ([57| in which a Langevin de- 
scription with local memory kernel (and thus, in the 
present case, also a white, Markovian noise as explained 
in the previous section) provides a valid description in 
the present model. The power law contribution p povi is 
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at most of relative order 7/M as compared to the expo- 
nential one in (|57|) for times t > tp with, in the limit 
where jt p ^> 1, 



7* 



where we used IjSoj) as well as Eq. 
pling g 2 for the damping rate 7. 
as the Langevin time. Demandin, 
on the power law contribution in 
more constraining condition (t F < 




(92) 

[86]) to trade the cou- 
We shall refer to t p 
; a similar constraint 
leads to a slightly 



t L p for M > 2m) 



hi 



n{/3M) 



n(2f3m) 



(93) 



Let us quote the relatively simpler expressions of the 
Langevin time (|92[) in the following cases. At low tem- 
perature, 1 <C (3m <C (3M, 



In 



,3/2 , , 

m\ M 



m 



(94) 



At intermediate temperature, j3m <C 1 <C 



7^ « In 



x 3 / 2 Jl/T 



7 



T 



(95) 



At high temperatures, (3m -C /3M <C 1, 

, 3/2 

7# « In (- 



(96) 



It is interesting to notice that the Langevin time is gov- 
erned by the ratio 7/m and not j/M. 

We close this section by mentioning that in the case 
M < 2m, the on-shell damping rate 7 = such that one 
has formally t pF — > 00 at one- loop. In this case, where 
dissipative processes tp — » \\ are kinematically forbiden, 
the pole contributions in Eqs. (|57|) - (|55|) are not damped 
and the power laws never dominates, hence t L p F — > 00. 
In fact, in this case damping arises at two-loop order. 
One has 7 cx <? 4 and the Langevin time is parametrically 
reduced by a In g 2 ~ In y/j contribution. 



V. CONSISTENCY CHECK 

We have shown an explicit example where the rele- 
vant self-energies exhibit a highly non local behavior, but 
where a local Langevin desciription still provides an accu- 
rate description in the time range It is instructive, 
for a deeper understanding of how this comes about as 
well as for discussing more intricate situations where an 
exact solution of the non-equilibrium equations of motion 
is not available, to analyze this issue directly in real time, 



at the level of the memory integrals in Eqs. (f27|) - (l28l) . in 
the spirit of our discussion in Subsec. Ill CI To this aim, 
we shall perform the following consistency check, com- 
ing back to the general discussion of Sec. Mil We plugg 
the solutions of the effective local description obtained 
in Subsec. IIII CI back in the original memory integrals of 
Eqs. ([27l) -(f28 |) and check whether the power law decay 
of the memory kernel ££(£) is sufficient for the validity 
of the local approximation (f3"Tj) . 

We begin with Eq. (f28|) or, equivalently, Eq. (|38|) for 
the spectral function p p (t). It is convenient to write the 
effective solution as 



Pf(t) = 



sin e„t 



.-7*1*1 = -ll m ( ( 

fir, L 



— ie p t--y p \t 



}, (97) 



which leads, for t > 0, to: 

pf (t) + Ippf (t) = Re {e-^+>)*} . (98) 

Expressing p p (t — r) in terms of p p (t) and dtp p (t) = p p {t), 
we get 

pf (* - T )= pf(t)Re [e^] - Pf® + W?® lm [e ~] , 

(99) 

with z = ie p + 7 p . Using this expression to evaluate the 
memory integral in Eq. (|63[) . we obtain 



drmt-T)pf(r) = I drmT)p; tt (t-r) 



„cff / 



where 



pf(t)+l P pf(t) 



a z p (t)= fdrW p {r)e^. 
Jo 



Ima-^(i), (100) 



(101) 



We are interested in the large time behavior of the 
integral (I101[) . We already see that for too large time the 
power law suppression due to the memory kernel is not 
enough to overcome the exponential e zt cx e 7p * and the 
integral diverges exponentially, signalling the breakdown 
of the effective description (|97l) . To evaluate the late time 
(e p t ^ 1) behavior of the integral (jlOll) more precisely, 
we introduce a time separation 77, such that l/e p <C rj <C 
1/7 P and rj <^ t, and write 

a;(t) « [ V dr S^(r)e"^ + [ dr X p p (r)e zr , (102) 

JO J-q 

where we neglected the exponential growth in the first 
term on the right-hand side. For t p r\ 3> 1, the upper 
bound contribution to the latter is suppressed by the os- 
cillating factor and one has 



dTE p p (T)e l ^ T 



dr^ p (r)e^ = ^(e p ). 



(103) 
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This is the required contribution for Eq. (|100[) to reduce 
to the desired local description, up to 0(j p /e p ) contribu- 
tions, see Eq. (|108[) below. 

To evaluate the second term on the right-hand side of 
(|102l) . we use the large time, power law behavior of the 
memory kernel, which we parametrize as 



a(pt) , 



(104) 



where a p is a dimensionful amplitude parameter, typ- 
ically proportional to some positive power of the rele- 
vant coupling constant, /i is a characteristic mass scale, 
a(x) ~ 1 is a bounded (e.g. periodic) function and v a 
given exponent [85| . We show in Appendix [B] that, for 

\zt\ — > CO, 

£*,, (T) ^£^£(-ff)-^ M . W 

where a^ n \x) — d n a(x)/dx n . Recalling the definitions 
(|32p and (|33"|) . we finally obtain, for e p t 3> 1, 



and 



Ima*{t) =ImE«(e p ) 



1 + 



1 + 



a p e 73 >* 



(106) 



(107) 



As long as the last terms in brackets in the above equa- 
tions can be neglected, the memory integral (|100|) has the 
desired form: 



dr££(t - r)pf(T) a te^t) + 2 7p/ >; tt (i) , (108) 



2 „eff / 



up to C(7 P /e p ) corrections, with 

7p = ' 



5e 2 p = Ret^(e p ) 



2e„ 



(109) 



The ansatz (JHTJ) is indeed a solution of Eq. ([351) witn CM]) 
and the local approximation is thus consistent within a 
finite time range governed by the terms in brackets in 
Eqs. (|106[) - (|107[) . For the local description - which in 
the present case gives rise to an exponentially decaying 
solution - to make sense, the latter should be larger than 
'Yp 1 which typically requires a weak coupling situation. 
The ratio a p j e p 5e p is parametrically of order one in the 
relevant coupling, whereas the ratio a p /j p e p may either 
be of order one in the case where damping is allowed at 
lowest order in perturbation theory, or be enhanced by 
inverse powers of the coupling if damping is only possi- 
ble at higher orders. Thus the most stringent restriction 
comes from Eq. (I107P and we check that the Langevin 
time, for which the local description makes sense, i.e. 



l P t L 




(110) 



agrees with the analysis of the previous sections, see 
Eq. (|57]). or Eq. flU) for the model discussed in Sec. [TV] 
A similar - although more lengthy - analysis can be 
performed for the -F-equation with the same conclusion, 
namely that, within the time range (|110p . 

/>oo 

/ dr W p (r)Ff(t-T, t')^Se 2 p Ff(t, t')+2 lp d t Ff (t, t'), 
Jo 

(111) 

up to corrections of relative order 0(7 P /e p ), where 
F° s (t,t') is given by Eq. (g7J) with F^(t£) given by flS5J) 
and Fp(t-t') replaced by F^ p (t-t r ) in §7$. This com- 
pletes our consistency check for the full non-equilibrium 
solution. 



VI. APPLICATION: DECOHERENCE, 
ENTROPY PRODUCTION AND 
THERM ALIZATION 

We end this paper with a simple application of the pre- 
vious considerations to the physics of quantum decoher- 
ence, entropy production and thermalization in quantum 
field theory in the context of the so-called incomplete de- 
scription picture, recently advocated in (58l-l65|. This is 
based on the observation that one's ability to measure 
the state of a quantum field is limited because of a re- 
stricted access to the infinite tower of n-point correlation 
functions: One typically has only access to low order cor- 
relators, most often to the subset of one- and two-point 
functions. The lack of knowledge of higher-order corre- 
lators may result, from the point of view of the observer, 
in effective loss of quantum purity and/or coherence and 
associated entropy production, even to effective thermal- 
ization |59l |60|. 

Let us consider the case where only the subset of inde- 
pendent equal-time two-point functions 



F p (t)= F p (t,t%, =t , 
R p (t) = d t F p (t,t%, =t , 
K p (t) = d t d t >F p (t,t%, =t 



(112) 



of a given field mode is measured at each time. From this 
knowledge one can reconstruct the least biased quantum 
state compatible with the measured correlators and infer 
its quantum properties. The corresponding density oper- 
ator is a Gaussian in the field operators, characterized by 
the intrinsic - canonically invariant - occupation number 

n p {t) + ^ = ^F p {t)K p {t)-Rl{t). (113) 

The latter measures the quantum purity of the Gaussian 
state (HI: 



P P (t) = 



1 



2n p (t) + 1 



(114) 



which is equal to its maximum value 1 (n p = 0) for a 
pure state. The occupation number (IL L3[) also measures 
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the phase space area covered by the Gaussian state in the 
Wigner representation and can be related to the Gaussian 
entropy [60j 

s p [t) = [n p (t] + l]\n[n p (t) + 1] - n p (t)\nn p (t) , (115) 

which measures the amount of missing information in the 
subset of measured correlators, that is the amount of in- 
formation stored in unmeasured higher-order correlators. 

Non-intrinsic - i.e. basis dependent - properties of the 
inferred Gaussian quantum state can be characterized by 
introducing another occupation number (note that < 
n p (t) < n p (t)) 



... 1 K p (t) + , 
np{t) + 2 = 



2 p F p (t) 



(116) 




For instance, the squeezing, or coherence parameter 



FIG. 1. The Gaussian purity (fTH)) (left) and entropy (TITS)) 
(right) as functions of x = 2 r y p t for the case of an initial 
vacuum-like - i.e. pure but with no long range coherence 
in the basis of coherent states - state: n p (0) = n p (0) = 0. 
The equilibrium occupation number, which determines the 
asymptotic values of both the purity and entropy (dashed 
lines), is n P q = 3. 



9 P (t) 



1- 



n p (t) + 1/2 
n p (t) + 1/2 



(117) 



measures the degree of quantum entanglement / coherence 
in the basis of semi-classical coherent states. A state 
with <7 p —¥ 1 exhibits non-trivial correlations between 
macroscopically distant semi-classical states and thus a 
high degree of quantum coherence (in this basis) [HI, [5t| . 

Now to the dynamics. We consider the simple situa- 
tion of Sec. IIIII where the field of interest is weakly cou- 
pled to a thermal bath with negligible backreaction. The 
exact solution of the nonequilibrium dynamics is given 
by Eqs. (|I7 j) -([i9 | . with ([53]). In the small width limit, 
the Langevin time (|57j) is large and the power law con- 
tribution in Eqs. (|55j) and ([5*5 |) - (j6"tJj) can be neglected, 
resulting in a Markovian dynamics as explained in Sub- 
sec. IIII CI Reshuffling the various terms (hence the bar 
on Fp below), one can write, for t, t' > 0, 



F p (t,t') =Fp(t-t') ri-e->[( t+t ')-|*- t 1] 
+i^ (M')e~ 7p(t+t ' ) , 



(118) 



where 



F°(t, t') =a+cos e p (t-t')+a p cos e p (t+t')+a° p sin e p (t+t'), 

(119) 

with, neglecting 0{^f p /e p ) corrections and setting Z p ss 1, 



F p (0) ± 



K p (0)\ »_R P (0) 



(120) 



In the Breit- Wigner - or Markov - approximation con- 
sidered here, the equal-time equilibrium two-point corre- 
lators are characterized by the on-shell equilibrium occu- 
pation number n p q — n(/3e p ): 



F^(0) 



1 



1 



(121) 



with i?^(0) = and K^(0) = e 2 p F°%0). It is an simple 
exercice to compute the various quantities (|113[) - (|117l) . 
One obtains 



n P (t) + \ = J[N?(t)] 2 + 2N e p «(t)N°(t) + [N0(t)}\ 

(122) 



where we defined 



N?(t) = (V + ^) (1 - e- 2lpt ) , (123) 
N$®= (^(O) + 0e- 2 ^ 4 (124) 



and 



JV°(i)= (n p (0) + ^je- 2 ^. (125) 

The occupation number n p (t) and the decoherence pa- 
rameter g p (t) have remarkably simple expressions: 



and 



i,(t) = (1 - e- 2 -"') + ft p (0)e- 2 >' 



n p (t) + 1/2 



(126) 



(127) 



Note that, as expected, an initially incoherent, or 
thermal-like state with g p (0) — or, equivalently n p (0) — 
n p (0), remains so: g p {t) = 0, or n p (t) — h p (t). A partic- 
ular case is n p (0) = n p (0) = 0, which corresponds to a 
pure, vacuum-like state. Preparing such a state at time 
t = 0, one observes an effective loss of quantum purity 
due to the incomplete knowledge of higher order corre- 
lators. Indeed, in that case, the Gaussian purity (|114[) 
reads 



P P (t) = T 



pcq 



(1 - P p cq )e- 2 ^* ' 



(128) 
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where P® q = l/(2n° q + 1). This loss of purity naturally 
corresponds to a loss of information and thus to a growth 
of entropy. Fig. [T] shows the purity (|114j) and the entropy 
(|115l) as a function of x — 2j p t for this case. 

Let us now consider the case of an initial pure state, 
n p (0) = 0, with a high degree of quantum coherence 
<7 P (0) — !• 1, i.e. n p (0) > 1, which corresponds to a 
squeezed vacuum-like state. Again because information 
about the state of the systems spreads towards higher 
order correlators, one observes apparent loss of quantum 
purity and associated entropy growth as well as apparent 
loss of quantum coherence. This is illustrated in Figs. [H 
and[3]for = 3 and n p (0) = 1, 3, 10. 

The case n p (0) = rip" 3 , which implies that n p {t) = np^ 
is particular. For instance, one gets 

(129) 

and thus the Gaussian purity and entropy approach their 
late time equilibrium value at a rate Ajp, twice the rate 
of the case n p (0) ^ rip* 1 . This is illustrated in the insert of 
Fig. [21 Also the coherence factor has a purely exponential 
decay law 

9 P (t) = S P (0K 2 ^ , (130) 

shown in Fig. [31 

For the case n p (0) > n,p q purity rapidly falls off and 
undershoots - entropy correspondingly overshoots - its 
equilibrium value before exponentially approaching it at 
a rate 2~/ p . The time i cxt at which an extremum is 
reached (minimum of P p or maximum of s p ) decreases 
when n p (0) is increased but remains bounded: 27 p i cxt > 
In 2. Thus the time scale of purity loss/entropy growth 
is always set by 7" 1 . Observe finally that cases with 
n p (0) > np q have faster loss of purity or growth of en- 
tropy than cases with n p (0) < rt° q , but they have slower 
loss of coherence. 



VII. SUMMARY AND CONCLUSION 

To summarize, we have seen that the dynamics of non- 
equilibrium quantum fields can always be formally writ- 
ten as an equivalent, non-local, non-linear Langevin pro- 
cess with self-consistently determined memory and noise 
kernels. A drastic simplification is achieved when the 
memory integral can be approximated by a local damping 
term. It is often assumed that this requires the existence 
of a finite, short memory time scale beyond which the 
memory kernel essentially vanishes. The typical situation 
in QFT is, however, that the latter exhibit a power law 
behavior at large times. Still, such a power law decay, if 
not enough to completely localize the memory integrals, 
is sufficient to make the dynamics essentially local in a 
finite time interval - before late-time power laws come to 
dominate - at least if the associated local damping rate 
is weak enough. 




FIG. 2. The Gaussian purity (TTH)) (left) and entropy (fTTSj) 
(right) as functions of x = 2-y p t for the case of an initial 
squeezed vacuum-like - i.e. pure and highly coherent - state: 
n p (0) = and n p (0) — 1, 3, 10 (top to bottom in the left panel 
and bottom to top in the right one), corresponding to initial 
coherence factors 1 — g p (0) — 5x 10 -2 , 10 -2 , 10 -3 respectively. 
The equilibrium occupation number, which determines the 
asymptotic values of both the purity and entropy (dashed 
lines), is n P q = 3. The insert shows the approach to the 
equilibrium value P p q in logarithmic scale. One sees that the 
relaxation rate is twice faster in the case n p (0) = n P q . The 
entropy, being one to one related to purity, exhibits a similar 
behavior. 



We have illustrated the above points in a simple ex- 
ample with non-equilibrium degrees of freedom interact- 
ing with a stationary thermal bath with negligible back- 
reaction, where an analytical solution of the QFT prob- 
lem is possible. We believe our analysis clarifies the valid- 
ity of local descriptions of non-equilibrium field dynam- 
ics, as well as of related Breit-Wigner motivated approx- 
imations, employed in a wide range of problems, see e.g. 
[3 GJ, H, [H H3|. Finally, we have provided a general 
method to work out the local limit directly at level of 
the equations of motion, which is useful when a complete 
solution is not known, for instance in non-stationary sit- 
uations, or when back-reaction is important - in which 
case the damping rate as well as the noise kernel may de- 
pend self-consistently on the non-equilibrium correlators. 

As an illustration, we have applied the above ideas to 
the description of quantum decoherence and entropy pro- 
duction in QFT in a simple situation with a stationary 
thermal environment along the lines of Refs [5{| [6(| • In- 
teresting further applications of the methods developed 
here include, for instance, the study of decoherence of 
primordial fluctuations in inflationary scenarios [58j], or 
the discussion of local damping descriptions in expanding 
geometries, which underly e.g. warm inflationary scenar- 
ios [25|, |86[ or the standard perturbative description of 
reheating in the early Universe [87j . 
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FIG. 3. The coherence parameter (I117P as a function of x = 
2y p t in linear (left) and logarithmic (right) scales for the same 
parameters as in Fig. f2J here n p (0) = 1, 3, 10 bottom to top. 
The case n p (0) = n° q leads to a purely exponential decay. 



Appendix A: Calculation of <r 



From Eqs. ([77 ]l -([78 ]l one easily obtains, after some 
manipulations, 



El 
~8tt 



dx\ 1 



2 m 



4m 2 / /3x 
— K~n 



fix 



(Al) 

Using n'(x) = —n(x)(n(x) + 1), one obtains the alterna- 
tive, simpler expressions 



9 

~4tt/3 
9 2 



+oc n((3mx) 

ax -j^=^= 

x 2 \/x 2 — 1 

l 



dx —= 

47T/3 Jo y/l 



X 



: n 



^ ] . (A2) 



X' \ x 

For instance, one obtains, at low temperature /3m ^> 1 



g 2 e -f)m 



inf3 



pra 



whereas in the high temperature limit j3m <C 1 

-,2 



(T 



9 



[1 + (/3m)] 



(A3) 



(A4) 



It is interesting to see how this constant arises in the 
standard calculation in frequency space, where one first 
computes Sp =0 (w) and deduce £^ =0 (w) by using the 
fluctuation-dissipation relation (|4"5]) . This seems to leave 
no space for the required S(u>) contribution in £^ =0 (u;). 
In fact, a careful evaluation of £p =0 (w) reveals a contri- 
bution 



(n(/3fc ) + i+n(/3(w-fc )) + ^ 



6(u) , (A5) 



where ko is a frequency variable to be integrated over. 
This contribution can be safely put to zero in T,p =Q (uj), 
owing to the property n(— x) + 1/2 = — (n(x) + 1/2). 
However, it leads to a non-vanishing contribution when 



multiplied by n{(3uj) + 1/2 to compute £p =0 (u;) as de- 
manded by the fluctuation dissipation relation. Indeed, 
using the identity 



1 

+ 2 



n(fi(u - ko)) 



n(/3w) 



n(/3(w - ko)) 



1 
2 
1 

" 4 ' 
(A6) 

+ 1 = 



which is easily demonstrated by using 2n(x) 
tanh~ (x/2) and making use of standard trigonometric 
relations, one obtains a non- vanishing term in £^ =0 (cj): 



- n(pko) 



6(u). 



(A7) 



One easily checks that this term is also obtained by a 
direct calculation of £p_ (cj), without using the fluctu- 
ation dissipation relation. Finally, one checks that the 
final coefficient of 5(uj) is precisely the constant <7oo ob- 
tained by the direct calculation in the time-momentum 
representation described above. 



Appendix B: Asymptotic behavior of (Xp(t) 



This section is devoted to the asymptotic behavior of 
J*dT~EP(T)e ZT in Eq. (fl~02|) as \zt\ -> oo. Integrating n- 
times by part reads: 

n ~ 1 (-~\) k 

dr^ p (r)e ZT = ]T {e zt s k (t) - e^s k ( V )} 

■n k =o 



+ (-l) n / dTS n { T )- 



(Bl) 



where s n {t) = ££(£). We recall that in the generic 
case considered here is a power law and therefore 

C°°. Thus for any finite interval [r),t], f dr s n (r)^- -+ 
as n — > oo. This allows us perform an infinite amount of 
integration by parts which gives 

(— 1 ) n 

dr££(r)e* T = £ ^ (e zt S „(t) - e*"s n fa)} . 

n>0 

(B2) 

Moreover, since Re z > 0, the contribution from the 
lower boundary is exponentially suppressed. Using the 
generic ansatz (|104l) and the fact that a is bounded, we 
get: 



where 



F{x,y) = Y J {-y) n ° {n) {x)V- + 0{l/x)] . (B4) 



n>0 
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The above expression simplifies when the function a(x) 
is periodic and thus satisfies a^ 2p ^{x) = (— l) p a{x) and ( f 

a (2 P +i)^ = (_i)P (T '( a; ). Splitting the sum into odd and f dl -Y,P( T )e ZT ~ — - (B6) 

even terms, one gets, for \y\ < 1, J v P (m*)" z 2 + y 2 

In the case where fi -C e p this further simplifies to 
T{x,y)= a{x) - y ^ x) [l + 0{l/x)\. (B5) z zt 

As \vl/z 2 \ = n 2 p /(e 2 p + 7p) < 1 we get 



hep-ph/0302210; hep- 
Conf.Proc. 739, 3 (2005), 



hep- 



J. Berges and J. Serreau, 
ph/0410330. 
J. Berges, AIP 
ph/0409233. 

J. Berges, S. Borsanyi, D. Sexty, and I.-O. Stamatescu, 
Phys.Rev.D 75, 045007 (2007), hep-lat/0609058. 
S. Borsanyi and M. Hindmarsh, Phys. Rev. D 79, 065010 
(2009), 0809.4711. 

C. Bodet, M. Kronenwett, B. Nowak, D. Sexty, and 
T. Gasenzer, (2011), 1101.0397. 

G. Aarts and A. Tranberg, Phys.Rev. D77, 123521 
(2008), 0712.1120. 

A. Tranberg, JHEP 0811, 037 (2008), 0806.3158. 

T. Prokopec, M. G. Schmidt, and S. Weinstock, Annals 

Phys. 314, 208 (2004), hep-ph/0312110. 

T. Prokopec, M. G. Schmidt, and S. Weinstock, Annals 

Phys. 314, 267 (2004), hep-ph/0406140. 

A. Anisimov, W. Buchmiiller, M. Drewes, and S. Men- 

dizabal, Phys.Rev.Lett. 104, 121102 (2010), 1001.3856. 

A. Anisimov, W. Buchmiiller, M. Drewes, and S. Men- 
dizabal, Annals Phys. 326, 1998 (2011), 1012.5821. 

M. Garny, A. Hohenegger, and A. Kartavtsev, (2010), 
1005.5385. 

B. Garbrecht and M. Herranen, Nucl.Phys. B861, 17 
(2012), 1112.5954. 

P. Arnold., G. D. Moore, and L. G. Yaffe, JHEP 0301, 
030 (2003), hep-ph/0209353. 

D. Boyanovsky, H. de Vega, and S.-Y. Wang, Phys.Rev. 
D61, 065006 (2000), hep-ph/9909369. 

Y. B. B. Ivanov, J. Knoll, and D. Voskresensky, 

Nucl.Phys. A672, 313 (2000), nucl-th/9905028. 

J.-P. Blaizot and E. Iancu, Phys.Rept. 359, 355 (2002), 

hep-ph/0101103. 

A. Anisimov, W. Buchmiiller, M. Drewes, and S. Men- 
dizabal, Annals Phys. 324, 1234 (2009), 0812.1934. 
M. Drewes, S. Mendizabal, and C. Weniger, (2012), 
1202.1301. 

M. Herranen, K. Kainulainen, and P. M. Rahkila, JHEP 
0905, 119 (2009), 0812.4029. 

A. A. Starobinsky and J. Yokoyama, Phys.Rev. D50, 
6357 (1994), astro-ph/9407016. 

T. Biro, Z. Xu, and C. Greiner, p. 226 (1998), arXiv:hep- 
ph/9809461. 

D. H. Rischke, Phys.Rev. C58, 2331 (1998), nucl- 
th/9806045. 

D. Bodeker, Phys.Lett. B426, 351 (1998), hep- 
ph/9801430. 

A. Berera and R. O. Ramos, Phys.Rev. D71, 023513 
(2005), hep-ph/0406339. 



[26] J. Yokoyama, Phys.Rev. D70, 103511 (2004), hep- 
ph/0406072. 

[27] J. Weenink and T. Prokopec, (2011), 1108.3994. 

[28] M. Nahrgang, S. Leupold, C. Herold, and M. Bleicher, 

Phys.Rev. C84, 024912 (2011), 1105.0622. 
[29] L. Kadanoff and G. Baym, Quantum Statistical Mechan- 
ics (Westview Press, 1994). 
[30] P. Danielewicz, Annals Phys. 152, 239 (1984). 
[31] E. Calzetta and B. Hu, Phys.Rev. D37, 2878 (1988). 
[32] M. Herranen, K. Kainulainen, and P. M. Rahkila, JHEP 

0809, 032 (2008), 0807.1435. 
[33] C. Fidler, M. Herranen, K. Kainulainen, and P. M. Rahk- 
ila, JHEP 1202, 065 (2012), 1108.2309. 
[34] S. Juchem, W. Cassing, and C. Greiner, Phys.Rev. D69, 

025006 (2004), hep-ph/0307353. 
[35] S. Juchem, W. Cassing, and C. Greiner, Nucl.Phys. 

A743, 92 (2004), nucl-th/0401046. 
[36] J. Berges and S. Borsanyi, Phys.Rev. D74, 045022 

(2006), hep-ph/0512155. 
[37] M. Lindner and M. M. Midler, Phys.Rev. D73, 125002 

(2006), hep-ph/0512147. 
[38] M. Lindner and M. M. Miiller, Phys.Rev. D77, 025027 

(2008), 0710.2917. 
[39] B. Garbrecht and M. Garny, Annals Phys. 327, 914 

(2012), 1108.3688. 
[40] R. Feynman and F. Vernon, Annals Phys. 24, 118 (1963). 
[41] M. Gleiser and R. O. Ramos, Phys.Rev. D50, 2441 

(1994), hep-ph/9311278. 
[42] C. Greiner and B. Miiller, Phys.Rev. D55, 1026 (1997), 

hep-th/9605048. 
[43] C. Greiner and S. Leupold, Annals Phys. 270, 328 

(1998), hep-ph/9802312. 
[44] D. Boyanovsky, K. Davey, and C. Ho, Phys.Rev. D71, 

023523 (2005), hep-ph/0411042. 
[45] N. Cassol-Seewald, R. Farias, E. Fraga, G. Krein, and 

R. O. Ramos, Physica A391, 4088 (2012), 0711.1866. 
[46] J. Rau and B. Miiller, Phys.Rept. 272, 1 (1996), nucl- 

th/9505009. 

[47] A. Berera, I. G. Moss, and R. O. Ramos, Phys.Rev. D76, 

083520 (2007), 0706.2793. 
[48] Z. Xu and C. Greiner, Phys.Rev. D62, 036012 (2000), 

hep-ph/9910562. 
[49] R. Farias, R. Ramos, and L. da Silva, Phys.Rev. E80, 

031143 (2009), 0909.2671. 
[50] R. Farias, R. Ramos, and L. da Silva, Comp. Phys. Comm. 

180, 574 (2009), 0909.2876. 
[51] E. Fraga, G. Krein, and L. Palhares, arXiv:0910.4369. 
[52] D. Boyanovsky, M. D'Attanasio, H. de Vega, R. Hol- 

man, and D.-S. Lee, Phys.Rev. D52, 6805 (1995), hep- 



15 



ph/9507414. 

[53] D. Boyanovsky, H. J. de Vega, R. Holman, S. P. Kumar, 
and R. D. Pisarski, Phys.Rev. D58, 125009 (1998), hep- 
ph/9802370. 

[54] D. Boyanovsky, H. de Vega, R. Holman, D. Lee, and 
A. Singh, Phys.Rev. D51, 4419 (1995), hep-ph/9408214. 

[55] L. Fonda, G. Ghirardi, and A. Rimini, Rept.Prog.Phys. 
41, 587 (1978). 

[56] A. Peres, Annals Phys. 129, 33 (1980). 

[57] H. Nakazato, M. Namiki, and S. Pascazio, 
Int.J.Mod.Phys. BIO, 247 (1996), quant-ph/9509016. 

[58] D. Campo and R. Parentani, Phys.Rev. D78, 065044 

(2008) , 0805.0548. 

[59] A. Giraud and J. Serreau, Phys.Rev.Lett. 104, 230405 

(2010) , 0910.2570. 

[60] J. F. Koksma, T. Prokopec, and M. G. Schmidt, 

Phys.Rev. D81, 065030 (2010), 0910.5733. 
[61] J. F. Koksma, T. Prokopec, and M. G. Schmidt, Annals 

Phys. 326, 1548 (2011), 1012.3701. 
[62] J. F. Koksma, T. Prokopec, and M. G. Schmidt, 

Phys.Rev. D83, 085011 (2011), 1102.4713. 
[63] T. Prokopec, M. G. Schmidt, and J. Weenink, (2012), 

1204.4124. 

[64] J. F. Koksma, T. Prokopec, and M. G. Schmidt, Annals 

Phys. 325, 1277 (2010), 1002.0749. 
[65] F. Gautier and J. Serreau, Phys.Rev. D83, 125004 

(2011) , 1104.0421. 

[66] J. S. Schwinger, J.Math.Phys. 2, 407 (1961). 

[67] L. Keldysh, Zh.Eksp.Teor.Fiz. 47, 1515 (1964). 

[68] J. Serreau, JHEP 0405, 078 (2004), hep-ph/0310051. 

[69] S. Borsanyi and U. Reinosa, Phys.Rev. D80, 125029 

(2009) , 0809.0496. 

[70] M. Garny and M. M. Miiller, Phys.Rev. D80, 085011 

(2009), 0904.3600. 
[71] For discussions of non-Gaussian sources, see [3(| H2, l7(il . 

[72] In general one should add the general solution of the 
homogeneous equation Gr*F = but the latter typically 
decays with time and has an infinite time to decay away 
in the present set-up where initial conditions are prepared 
at t — > — oo. 

[73] Initial conditions for the spectral function p are fixed by 
equal-time commutation relations. 

[74] This (exact) equivalence has been repeatedly noticed in 
the literature for the case where the field of interest is 
linearly coupled to an environment on which it does not 



backreact, see e.g. [T^, |40| . [44| . In more general cases, 
the influence functional approach of Feynman and Ver- 
non allows one to derive approximate Langevin equations 
[23l. |40|. l47l] . To our knowledge, the - formal but exact - 
equivalence described in the present paper has not been 
pointed out before. 

[75] G. D. Moore, (2002), hep-ph/0211281. 

[76] This is somewhat similar to so-called harmonic, or one- 
frequency ansatz employed in Refs. [23|, 52, IHf. Here, 
the convolution with the memory kernel E p controls the 
convergence of the integral and gives better control of the 
approximation. 

[77] It is also useful to recall that G>(-w) = [Gf(ui)]*, 
Aj(-w) = -ft,(w) = \p P (uj)r, F^(-oj) = F^(oj) = 
[F p q (uj)]* and similarly for the various self-energy com- 
ponents. 

[78] M. Drewes, (2010), 1012.5380. 

[79] There is a similar expression near the symmetric pole 
to — > — e p - recall that p p {— lo) — —p~ p (u>). 

[80] All these quantities may in principle depend on the mo- 
mentum p. We omit the p index for simplicity. 

[81] The solution obtained here generalizes that of Ref. [3!|, 
which corresponds to the particular nonequilibrium con- 
ditions a~ — a p = 0. 

[82] We choose this ansatz for convenience although the poles 
are slightly displaced by a negligible amount ~ 7p/ e p 
along the real axis. In particular, it exactly reproduces 
Eq. (|50jl - with no j p correction - when u> — > e p . 

[83] We emphasize that this constant only appears for the 
mode p = 0. Also, it is a particular feature of the present 
model with cubic interaction vertex. For instance, it is 
easy to check that there is no such constant term for the 
lowest non-trivial contribution to damping (a two-loop 
setting-sun diagram) for a theory with quartic interaction 
vertex. 

[84] We assume the small width limit 7/M <C 1 throughout 
and neglect the mass correction ~ ReEij(M) as well as 
the field renormalization (|51|l for simplicity. 

[85] All these quantities may in principle depend on the mo- 
mentum p. We omit the p index for simplicity. 

[86] G. Aarts and A. Tranberg, Phys.Lett. B650, 65 (2007), 
hep-ph/0701205. 

[87] R. Allahverdi, R. Brandenberger, F.-Y. Cyr-Racine, and 
A. Mazumdar, Ann.Rev.Nucl.Part.Sci. 60, 27 (2010), 
1001.2600. 

[88] J. Berges and J. Cox, Phys.Lett. B517, 369 (2001), hep- 
ph/0006160. 



